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Abstract 

A systematic study of deformations of four-dimensional Einsteinian space-times 
embedded in a pseudo-Euclidean space E'^ of higher dimension is presented. In- 
finitesimal deformations, seen as vector fields in i?^, can be divided in two parts, 
tangent to the embedded hypersurface and orthogonal to it; only the second ones 
are relevant, the tangent ones being equivalent to coordinate transformations in 
the embedded manifold. 

The geometrical quantities can be then expressed in terms of embedding func- 
tions and their infinitesimal deformations z^ = z^ + ev^. The de- 
formations are called Einsteinian if they keep Einstein equations satisfied up to a 
given order in e. The system so obtained is then analyzed in particular in the case 
of the Schwarzschild metric taken as the starting point, and some solutions of the 
first-order deformation of Einstein's equations are found. 

We discuss also second and third order deformations leading to wave-like so- 
lutions and to the departure from spherical symmetry towards an axial one (the 
approximate Kerr solution) 
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1 Introduction 



The study of the two-body problem in General Relativity must include possible emission 
of gravitational waves. In this respect General Relativity runs into difficulties akin to 
those of classical electromagnetism, where the full treatment of charged particle motion 
plus the radiative field can be performed only via approximation techniques. Great 
progress was achieved since the first papers by LI. Bel and N. Deruelle, [1], followed 
by the throughout analysis of post-Newtonian and post-post-Newtonian approximations 
including gravitational radiation. The most important contributions have been made by 
Th. Damour, L. Blanchet, G. Shaeffer and many others, [2], [3], [1]. 

Their approach was based on a particular choice of the initial approximation. In 
physical situations where the gravitational field is not too strong and when the relative 
velocities of bodies under consideration are very low as compared with the speed of light 
c, the Newtonian theory provides us with exact (although not explicit) solutions. One 
takes one of these solutions as a starting point, adding successive corrections resulting 
from the inclusion of the relativistic effects: the finite propagation of gravitational field, 
the curvature of space, and so forth. The corrections are consequently made to the 
trajectory, to the law of motion, and to the gravitation potential (identified with the 
corrections to the space-time metric tensor). 

In a series of papers published a few years ago [B], [T|, [S], [H]) an alternative method of 
determining relativistic motion of test particles in a spherically symmetric gravitational 
field has been proposed, without need to use the Newtonian limit of General Relativity. 
The idea was to take as a starting point the very special explicit solution of geodesic 
equation in Schwartzschild metric background: a test particle moving with constant 
speed along the circle around a spherically symmetric mass. Let us denote this particular 
trajectory by (s), with s denoting the proper spacetime length of the curve. Then one 
can consider a deviation from this worldline, (s), which can be expanded in a power 
series of some small parameter e: 



The geodesic equation in Schwarzschild background can be also expanded in a series of 
equations whose solutions in terms of unknown deviations ra^ (s), 6^ (s), etc., will provide 
us with successive approximations to the exact solution, up to any order required. The 
parameter e is roughly proportional to the eccentricity of the new trajectory. Already 
the first approximation, linear in e, predicts the perihelion advance for orbits with very 
small eccentricity. 

The shortcoming of this method, which was entirely focused on the trajectories, was 
the total lack of any variation of the gravitational field, i.e. the Schwarzschild metric 
which was maintained invariable for all orders of geodesic deviation. This fact reduced 
the validity of the method only to the case of test particles with mass m negligible 
when compared with the mass M of the central body appearing in the Schwarzschild 
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background metric. More precisely, the successive approximations of planet's trajectory 
remain valid as long as the dimensionless parameter (m/M) can be considered as negli- 
gibly small, i.e. (m/M) << 1. When the mass m is no more a negligible quantity, its 
presence must inevitably alter the geometry of the initial Schwarzschild metric, and its 
influence can be therefore represented by a power series in the small parameter (m/M). 

In this article we shall describe the departure from the initial Schwarzschild metric in 
terms of the embedding functions. Embeddings of the exterior Schwarzschild geometry 
in pseudo- Euclidean fiat spaces are known since a long time ([I3], [H], [I2]), and once 
such an embedding is given, all intrinsic geometric quantities of the embedded manifold 
can be expressed in terms of derivatives of the embedding functions which depend on 
four "internal" parameters which are the space-time coordinates. A general analysis 
of deformations of the embedded Einstein spaces was given in [10]; nevertheless, only 
the theoretical setup was considered, without any concrete solution describing Ricci-flat 
deformations of known exact Einstein spaces, and in the first place, Schwarzschild and 
Kerr metrics. 

The present article is intended to explore not only the first linear approximation, but 
also the effects of second and third order in the expansion of deformations in powers of 
small parameter e, including the corrections describing gravitational waves. The depar- 
ture from spherical symmetry (Schwarzschild's metric) towards axial symmetry (Kerr's 
metric) as Einsteinian (Ricci-fiat) deformation of the corresponding embedded manifold 
is also discussed. 

2 Isometric embeddings and their properties 

2.1 The embedding functions and the induced metric 

Consider the embedding of a four-dimensional Riemannian space parametrized by local 
coordinates (denoted by x^, /i, = 0, 1, 2, 3 as usual) in a pseudo-Euclidean space of 
dimension A^. The dimension N, yet unspecified, depends on the topology of the Rieman- 
nian space under consideration, and may be quite high, as acknowledged in [12] . Locally, 
any n-dimensional Riemannian manifold can be embedded in a (pseudo)-Euclidean space 
of dimension N = n{n + l)/2. Here we are interested in global embeddings, which 
may require a relatively low dimension of the "host" space if the Riemannian space to 
be embedded possesses some particular symmetry. For example, the de Sitter space 
can be embedded globally in a five-dimensional pseudo-Euclidean space with signature 

(H ), and both exterior and interior Schwarzschild solutions can be embedded 

globally in a six-dimensional E'^ with signatures {+-\ ) or (H ). Consider 

a global embedding of a Riemannian space V4 given by the following set of embedding 
functions z^: 



z 



A 




A 



with 



1,2, ...AT, 
0,1,2,3. 



(2) 
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The metric tensor of V4 is the induced metric defined as 



9^,u = VAB d^z^ (3) 

o 

The inverse metric tensor 9 cannot be obtained directly from the embedding func- 
tions, but should be computed from the covariant components as their inverse matrix. 
From now on we use the superscript notation in order to make difference between the 

o 

"basic" induced metric 9 which will be considered as a background, and its infinitesimal 
deformations expanded in terms of an infinitesimal parameter e as follows: 



^M^^ = V + ^V + ^^V + -- (4) 

induced by the following deformation of the initial embedding functions: 

z^ix'')^ z^{x'')+ev^{x>')+e^l^{x>') + ... (5) 

When seen from the ambient pseudo-Euclidean space, the new embedded manifold 
is the result of an infinitesimal deformation of the initial manifold V4 induced by a 
vector field in E^^y It is quite obvious that on the embedded manifold such a field 
can be decomposed into its normal part (in the sense of the pseudo-Euclidean metric) 
and a part tangent to V4. This last part induces an internal diffeomorphism of V4 and 
can be always implemented as a local coordinate transformation. Such deformations do 
not have any physical meaning, but it is not always necessary to consider exclusively 
the deformations orthogonal to the embedded V4; sometimes a deformation having non- 
vanishing both parallel and orthogonal parts can have less non-zero components in the 
ambient space Efi s than its part orthogonal to the embedded V4 manifold. 



2.2 Expressions for connection and curvature 

Our first aim is to express all important geometrical quantities e.g. the connection 
coefficients and the curvature tensor, in terms of embedding functions z^ and their partial 
derivatives. Let us start with Christoffel connection 

r = ^ ^ (^dj up + dj - d^g . (6) 

o 

/^From the definition oi 9 (13]) we have the expression for its partial derivatives: 

dj,. = VAB (^dl z^ z^ + d, z^ dl, zA^ . (7) 

When substituted into the definition ([6]) it gives 

vl. = ^AB>d,z^dl,z^ (8) 
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Consider now the second covariant derivative of z^: 



(9) 



Therefore, we have 



5>'' = v^v.^^ + ri^^AZ^, (10) 

and of course, z^ = z^ . Substituting ( fTOl) into (IE]) we obtain the following identity: 



9 V, z^ V^V, z"" + 7]AB z"" Vp z^. 



i- flU — VAB ^ Vp^ V^Vjy 

But in the last term we note that 

9 11 AB V, z^ Vpz"" = 9^'^hp = S^, 

reducing (fTTj) to 

f J. = VAB 9 ^"Vp z"^ Vp V, 2;^ + f 
which shows clearly that 

VAB^pZ'^V^V.z'' =0 



X 



fill 



(12) 



(13) 



(14) 



which may be considered as an alternative (although implicit) definition of Christoffel 
symbols, and could be also derived as a direct consequence of the fact that gxp = 0. 

Using this result, let us form the following combination of covariant derivatives which 
vanishes identically: 



VAB 



Vp( Vp V.V. z"") - V.(Vp z"" VpV, z 







Applying the derivation and using the Leibniz rule we get: 



VAB 



V^Vp V.V. - V. Vp z^ Vp V. z"") 



+ 



+ r]AB^pZ^ 



Vp V,V, 2;^ - V, VpV, z' 



(15) 



Recalling that 



VpV. -V.Vp 



Vp2;^ = i?p/pVA z"", 



(16) 
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so that we can write 



o 

R 



nv\p — —Vab 



(17) 



which is the well known Gauss-Codazzi equation. 

The definition of the Riemann tensor by means of derivatives of the embedding func- 
tions given by formula ( ITTI) looks very compact, but is in fact highly non-linear and quite 
complicated. This is so because it contains many Christoffel symbols involved in the 

o 

second covariant derivatives, which contain in turn the contravariant metric tensor 9^"^. 
The components of the contravariant metric tensor are obtained as rational expressions 
in third and fourth powers of z^. Nevertheless, the most important point here is that 
the Riemann tensor depends only on first and second derivatives of embedding functions, 
so that the Einstein equations expressed in terms of the embedding functions will lead 
to second-order partial differential equations. 

The expressions derived in this section will be very useful in the development of a 
power series expansion of infinitesimally deformed embedding. 



3 Infinitesimal deformations of embeddings 
3.1 General setting 

o 

Let us consider an isometric embedding of an Einsteinian manifold Va in a pseudo- 
Euclidean space Ep with signature {p+, q—), with p + q = N: 



Va 



A, 5,... = 1,2,. ..AT, 
P'? with ^\ o o 

z^(x^) = 0,1,2,3. 



Consider now an infinitesimal deformation of the embedding defined by a converging 
series of terms proportional to the consecutive powers of a small parameter e. The 
deformed embedding defines an Einsteinian space V4; 



(x^) ^ 5^ (x^) = 2^ (x'^) + £ t; ^ (x^) + e' w"" (x'^) + . . . 
The induced metric on V4 can also be developed in series of powers of e: 



Vab 



d^z^d.z^ +e( d^z'^d.v^ + d^v'^d.z^ ) + 



(19) 
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Among all possible infinitesimal deformations of the embedding functions (x^) + 



ev 



there is a large class of functions 



which will not alter the intrinsic 



geometry of the embedded manifold. Infinitesimal translations v"^ = Const, obviously 

O /ID 

do not change the internal metric 9 = f]AB dfiZ dvZ ■ Also the generalized Lorentz 
transformations of the pseudo- Euclidean space E^^-^ keep the internal metric unchanged. 
Indeed, if we set 



5^ 



A R 



(20) 



with A'^g constant matrix. Then the first-order deformed metric is: 



9iiu = 9 tj,u + £9 + ...^ 

= TlAB 



d^z^d^z^ + e (^d^z^A^d^z"" + k'^cd.z^d^z'' 



= 9^,u + e 



d^z'^d.z^ 



Then the first-order correction vanishes if the matrices A^ satisfy the identity 

VAB + VCB^^A = 

which defines the infinitesimal rigid rotations (Lorentz transformations) of the pseudo- 
Euchdean space Ej^^y 

The geometric character of our approach enables us to eliminate unphysical degrees of 
freedom using simple geometrical arguments. Remember that in the traditional approach 
leading to linearized equations for gravitational fields the starting point is the following 
development of the metric tensor: 



(21) 



thus introducing ten components of hf^u as dynamical fields. We know however that most 
of them do not represent real dynamical degrees of freedom due to the gauge invariance. 
The metric tensor itself does not correspond to any directly measurable quantity. In 
fact, its components may be changed by a gauge transformation without changing the 
components of the Riemann tensor which is the source of measurable gravitational effects. 
In particular, the gauge transformation 



(22) 



does not alter the Riemann tensor so that both gi^,^ and Qfj^v describe the same gravitational 
field. 
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The arbitrary vector field generating gauge transformation fl22|) represents four 
degrees of freedom which are redundant in g^^,; this is why in the hnearized Einstein 
equations one may impose four gauge conditions e.g. 

h"" = 0. (23) 

The unphysical degrees of freedom can be easily eliminated from the embedding de- 
formation functions f'^(x^) if we note that any vector field in the embedding space 
that is tangent to the embedded Riemannian space V4 describes nothing else but a dif- 
feomorphism of V4, in other words a coordinate change, which has no influence on any 
physical or geometrical quantities. 

Vector fields tangent to the four- dimensional embedded manifold V4 can be decom- 
posed along four arbitrarily chosen independent smooth vector fields in tangent to 
V4. On the other hand, vector fields transversal to the embedded hypersurface V4 must 
satisfy the following obvious orthogonality conditions: 

VAB = r]AB = 0. (24) 

For any value of A the four partial derivatives (let us remind that z"^ = d^^ z"^) span a 
basis of four vector fields in E^ tangent to the sub manifold V4; therefore any vector 
satisfying the orthogonality condition ([21]) is transversal to V4 (as seen in E^). 

The orthogonality condition (12^ imposes four independent equations, which reduce 
the number of independent deformation functions f to iV — 4. This means that gen- 
eral non-redundant deformations can be decomposed along N — 4 independent fields 
A; = 1,2,..., AT -4: 

Af-4 

.-^(a;^) = 5^t;'=(x^))X(t)(x^). (25) 

k=l 

The basic fields X^^ (x^) can be chosen at will provided they induce a non-singular global 
vector field on V4, while the relevant degrees of freedom are contained in — 4 functions 
f'^(x^). To take an example, the de Sitter space can be globally embedded in a five- 
dimensional pseudo- Euclidean space 4 with signature (H ); therefore its global 

deformations can be described by a single function v{x'^) (see p!0]). 

One may ask the following question: if the deformation destroys the initial symme- 
try of the embedded manifold so that the deformed manifold cannot be embedded in 
the initially sufficient A^-dimensional pseudo- Euclidean space but needs a fiat embedding 
space of higher dimension ? It is known that global embeddings of the Kerr metric 
need more than six fiat dimensions sufficient for the embedding of the exterior (or in- 
terior) Schwarzschild solution (see [19], [22], [23], [21]), although Schwarzschild's metric 
can be obtained from Kerr's metric as a limit when the Kerr parameter a (the angular 
momentum) tends to zero. 

The answer is that as long as we investigate only the first-order corrections to geom- 
etry, we should not worry about this issue for the following two reasons: first, when a 
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global embedding is given, its infinitesimal deformations cannot lead to a global mod- 
ification of the embedding; second, if the bigger embedding space was introduced, say 
£jN+m^ it would contain the initial embedding space as its linear subspace, so that 

^N+m eE"", 

and its pseudo-Euclidean metric could be represented as a blockwise reducible matrix 



(26) 



^Vab 
VijJ ' 

with A,B,... = 1,2,...N, = 1,2, ..m, a, (3,... = 1,2, N -\- m. Accordingly, any 

deformation of the initial embedding can be decomposed in two parts, one contained in 
the initial embedding space E^ and another one in the complementary subspace E^: 

v"^ ^[v^,v'^\. (27) 

But the initial embedding functions had their components entirely in the first subspace 
E^ , = [z-^, 0], therefore the deformed embedding functions can be written as 



z"" = [z"^ + ev^, ev^], 
so that the induced metric of the deformed embedding will be 

o 1 „2 

Qnu = 9 ^,u + £9 + £ 9 + - 

= r]AB (d^z''d,z''\ + e tjab (d^z^d.v'' + Oyd.z'' 



{2i 



(29) 



VAB ( d^v'^d.v'' + d^z'^d.w'' + di^w^'d^z'' ] + r]ij ( v'd, 



.3 



B 



^Prom this one can see that the deformations towards the extra dimensions do not con- 
tribute to the first-order corrections of any geometrical quantities obtained from the de- 
formed embedding functions. This is why we shall not consider such deformations while 
investigating at first only the terms linear in the infinitesimal parameter e. Our principal 
aim now is to establish the explicit form of connection and curvature components in- 
duced on the infinitesimally deformed embedding V4. To this end we must calculate the 
approximate expression of the contravariant metric tensor g'^'^. If the covariant metric is 
decomposed as 

^M'^ = V + ^V + ^^V + -" (30) 

then we have the following formulae defining the corresponding decomposition of g'^'': 
o 1 



9' 





2 






9 ^'9'"' 9 pa 





o o 0,11 



(31) 
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3.2 The first order corrections to Einstein equations 

In what follows we shall keep only the first order terms linear in e. 

Let us start by computing the first (linear) correction to the components of the 
Christoffel connection, which develops in Taylor series as 



r 



pv 



° A 1 A 2 , 

-T^ -x- eV +r + 

^ pu ' ^ ^ pu ' ^ ' 



(32) 



then by definition we have: 



pu 



1 o \ / 1 1 



1 

up pp dp9 pu 



1 o 

2 



1 



1 



Y^^^i dp9 up + 0^9 pp - dp9 



1 



- 9^' { dp9 up + d,9 pp - dp9 pu] - 9^" 9'''9 .J d,9 uo + d,9 



o 



o 



CTK 1 '-'pJ Up I '-'UJ pp 



pu 



do9 



(33) 



p^ pu 



One easily checks that 



1 o - / 1 1 

^9''{dp9up + d,9 



pp 



do9 



p^ pu 



Vab9 



Xp 



dpz^dl^v^' + dpv^dlz^ 



(34) 



while the first term after some algebra gives 



1 1 



9^^{ dp9 up + d^9 pp - dp9 pu 



-Vab9 



Xp 



dpZ^T%d^v'' + dpv''T%d,z'' 



(35) 



Combining together ( IMI) . ( IMl) and ( l35l) we find the final expression 



pu 



r]AB9 



Xp 



Vpz^ Vu V, + Vpv^ VpVu z^ 



(36) 



This expression has a tensorial character as it should be, because by definition both 
quantities 



r and r 



A 
pu 



transform as connection coefficients, therefore their difference must transform as a tensor, 

and this is true for any term of the development into series of powers of e. 

1 

The coefficients Tpy will be useful for the derivation of geodesic equations in the 
deformed space-time, but they are not necessary for the computation of the first-order 
deformation of the Riemann tensor, which can be determined as follows. 

Let us develop second covariant derivatives of the deformed embedding functions z^ 
yields: 



ypV,z^ + eVpVuV^ + 0{i' 



VpV, z'^ + e 



(37) 
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The Riemann tensor induced on the deformed embedding is defined by the same formula 
as in the previous section f|T7|) : 



RiyfiXp — Vab 



V^Va V, Vp 5^ - V,Va 5^ VpVp z 



A ■ 



{31 



Note that in order to calculate the components of the Riemann tensor induced on the 
deformed manifold V4 we use not only the deformed embedding functions z^, but also 



the "deformed" covariant derivations V^. 



Now, when we insert the expressions like ( 1371) into the definition of Riemann tensor 
components ( l38l) . we shall encounter, besides the zeroth-order initial Riemann tensor 

RfiuXp and the second-order corrections proportional to e^, just two types of terms linear 
in e: 



eriAsV^Nxz'^V.Vpv'' 



and 



.B 



up 



(39) 



The terms of the second type vanish by virtue of the identity ( IT^ : therefore the first-order 
correction to the components of Riemann tensor can be written as follows: 



R 



Vab 



Va z"" VuVo v"" + V„ Va V^Va z 



U V p 



pVA' 



u y p ■ 



(40) 



To establish the form of linear correction to Einstein's equations we need to know the 
components of the first-order correction to the Ricci tensor and the Riemann scalar. 
These quantities are readily computed as follows: 



R pp — 9 R pvxp puXpi 
Consequently, the first-order correction to the Riemann scalar is: 



(41) 



1 o 1 1 o 

R = 9'"' Rpu + 9'"' R 



(42) 



Finally, The first-order correction to the Einstein tensor, i.e. the left-hand side of Ein- 
stein's equations is: 



1 

G 



pu 



1 1 1 o 1 o 1 

R pv '^9 pu R ~ 9 pi/ R 

1 

R 



(43) 



1 o o - 1 



1 o 1,0 

Ap - 



pv ^2^9 9 R Xp r, 9 pu 9 Rxp r, 9 pu R- 



1 1 o 

2 
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In what follows, we shall always suppose that the initial Riemannian manifold is a 
solution of Einstein's equations, i.e. an Einstein space which is Ricci-flat and consequently 
has zero scalar curvature, too. Therefore the linear correction (of the first order in small 
parameter varepsilon ) to the Einstein tensor will reduce to: 

R^^u-Y^^.u9^'Rxp (44) 

In the absence of any extra gravitating matter (besides the matter generating the basic 
solution, e.g. the central spherical body for Schwarzschild's solution) the equations to 
solve can be written in form of a matrix acting on the first-order correction to the Ricci 
tensor: 

^^i^'u-lh^a'^^ Rxp = (45) 

But this amounts to the Ricci flatness up to the first order, because the operator acting 
on the right on the Ricci tensor in (jlSD is non-singular; in fact, it is its own inverse: 

StK-lLJ''^(^S;6:-^9,J^^^=6iS^^ (46) 

iFiom this we infer that in an Einsteinian background the first-order correction in vacuo 
should satisfy the equation 

Rxp = 0. (47) 
This may be written, developing ( l42l) . as: 

= Uj\V,V,v^ + Vj^Vy = (48) 

with 

Uj\ ^ vab9 {SlSJV.Vy + SiS:V,V,z^ - SiSJV.V.z^ - Sl52V.V,z^) (49) 

- 4v«^^V^V^z^V,Vp^^ - diVpz^'V^V^z'^V^Vpz''^ (50) 

In the case when the energy-momentum tensor is present (supposing however that it 
describes the influence of matter weak enough in order to keep the basic solution un- 
changed), one must use the full Einstein's tensor on the right-hand side. The first correc- 
tion, linear in e, reduces then to only two terms due to the fact that the initial solution 



is an Einstein space in vacuo so that Rxp = and R: 



1 1 o ° Xp'^ 

R fii/ '^9 i_iu 9 R Xp 



\ loo. 



R fiu — ~T^^, (51) 
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and this in turn, due to the idempotent property fH6l) . can be written equivalently as 

SttG 



1 

R 



Xp 



T: 



Xp 



(52) 



which may prove to be more practical for further calculations especially when the energy- 
momentum tensor has a particularly simple form. 

3.3 Second order deformations 

Expanding the deformed embedding functions into the power series of small parameter 
e as earlier seen: 



(53) 



the metric tensor g^^^ of the deformed embedding is given by the formula ([T9|) . 
The expansion of connection coefficients was already given in eq. ( l32l) . 
The second order terms in the expansion of the Riemann tensor are given by the 
following formula: 



2 

Rpuap 



Vab 



+ V.V.w^ - r,,V«i;^ - r,,V«2;^ V^VpZ^ - (/i ^ z/ 



(54) 



The terms proportional to V are zero according to (fT^ . With a little algebra we find: 



2 

Rpuap 



= VAB [VuV^z^V^VpW^ + V.V^w^V^VpZ^ - 



o l'^ 1-^ 



1 1 



The r r terms contain only z and v functions. 

To write the second-order correction to Einstein equations we need the correction to 
the Ricci tensor: 



o 2 



1 1 



2 o 



Rup 9 ^ Rp,uap ~\~ 9 ^ Rfiuap ~\~ 9 ^ Rfiuap 
2 2 

with w functions contained only in Rp^ap and in 5'^°'. 



(55) 
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Looking for vacuum solution we must develop this equation. For the sake of simplicity 
we note only that the operators acting on the derivative of w"^ are the same given in fHOl) 
and ( l50l) . The know functions and v"^ serve now as the right-hand side of the equations 
determining the w"^ functions: 



2 



i?.p = ^ Uj\V^V,w^ + Vj^V^w^ = B,,{v^, z^) 
with Bypi^v"^ , z"^) combination of derivative of v"^ and z"^ functions. 

4 Approximate solutions of Einstein equations 
4.1 Flat background space-time 

In a Minkowskian space-time M4 parameterized by cartesian coordinates = [ct, x, y, z] 
all connection coefficients identically vanish, as well as the components of the Riemann 
and Ricci tensors. The flat Minkowskian space can be embedded as a hyperplane in any 
pseudo- Euclidean space with more than four dimensions and signature (1+, (A^ — 1)— ). 
Let us choose the simplest case of embedding in five dimensions: 



M4 El 



5 



,4 



with the first four components denoting a Minkowskian space-time vector in cartesian 
coordinates: 

= ct, z^ =x, z^ = y, z^ = z, z^ = 0, (56) 

the last cartesian coordinate considered as an extra dimension of Ef^ orthogonal to the 
M4 hyperplane. All covariant derivatives in fl38|) can be replaced by partial derivatives, 
and all second derivatives of linear embedding functions are identically zero. Therefore 
in order to investigate non trivial deformations of the Minkowskian space embedded as 
a hyperplane we must go the second order in e. This leads to the following equation 
resulting from the requirement of vanishing Ricci tensor: 



2 oAi^ 
Rfip = =^ 9 TjAB 



V. Va V. Vp v"" - V. Va V. V„ v"" 







We shall not consider infinitesimal deformations of the first four coordinates because 
they coincide with coordinate transformations in V4; therefore the only non vanishing 
component of is the remaining fifth coordinate deformation, expanded in a series of 
powers of e: 

/ = ev^x^") + + + ... 
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In order to keep the Einstein equations satisfied after deformation up to tlie second order 
terms, we must have 



R. 



9 



Xu 



(57) 



Any function of hnear combination of cartesian coordinates is an obvious solution of Eq. 
fl571). Indeed, if we set: 



(5^ 



inserting the derivatives of v{k^x^) into fl57|) results in the following simple equation 



9 



Xu 



knkx kukp V kyk\ knkp v 



k^k" V 



i2 



kpkp kpk^ 



0, 



(59) 



But in fact, this deformation does not have any physical meaning, because the Riemann 
tensor, which is the only observable quantity, identically vanishes: 



2 

Rfiuxp 



kp, k\ ky kp ku k\ /c^ kp 







(60) 



The vanishing of the Riemann tensor is not surprising, because the deformation consid- 
ered looks like a deformation of a plane into a cylinder, which does not alter its intrinsic 
fiat geometry. 

The fact that there are no wave-like solutions at the first order of deformation of 
Minkowskian spacetime suggests that the same situation will prevail when we shall in- 
vestigate other Einsteinian manifolds embedded in a pseudo-Euclidan fiat space, e.g. the 
Schwarzschild solution. If the contrary was true, one could keep the wave-like propa- 
gating deformations also in the fiat limit, which would contradict the absence of such 
solutions among the first-order deformations of the Minkowskian space-time. 

This means that the only hope to produce contributions to the Riemann tensor be- 
having like a propagating gravitational field, i.e. the gravitational waves, is to consider 
the third (and higher) order deformations of embedded Einsteinian manifolds. The third 
order variation for the Riemann tensor in the case of deformations of all orders orthogonal 
to the embedded manifold reduces to the following expression: 



3 

Rp,uap 



VAB {VpVpV^V^VaW^ + VpVpW^VuVaV^- 
VMpV^VpV^w'' - V.VpW^VpV^v'') 



(61) 



The linear contribution coming from the expressions containing third-order devia- 
tion linearly does vanish because the derivatives of the corresponding coordinate are 
identically zero. 

A wave-like behavior of the Riemann tensor can be produced if we assume that w 
depends on variables orthogonal to the worldlines parallel to the vector k. For the sake of 
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simplicity, let us start with the first order deformation in the direction of fifth coordinate, 
i.e. orthogonal to the embedded Minkowskian hyperplane M4 as a plane wave propagating 
along the z-axis: 

According to our general analysis, by virtue of (jSO]), this deformation does not con- 
tribute to the Riemann tensor, which remains zero even at the second order. Now let us 
add up the second order deformation depending on the variables x and y only: 



= ee*("*-^") + e'ti;(x,i/) (62) 

The only contribution to the third order correction to the Riemann tensor has the form 
given by the formula ( 16T1) in which the covariant derivatives can be replaced by partial 
derivatives given that all Christoffel symbols vanish in cartesian coordinates. The func- 
tion w{x,y) must have some non vanishing second order derivatives; let us make the 
simplest choice and set w{x,y) = Bxy, with B = Const having the dimension cm~^. 

Then the only non vanishing second derivative is d^yW = B. Taking into account the 
form of (16T]) . the only non vanishing components are: 



Rt^xyp = {dly vdlpW + dly w dip v- 8% v dl^w - 8% w d^^ v) (63) 

and all other components obtained from this one by permutations of indexes allowed by 

3 

the well known symmetries of Riemann's tensor, like e.g. Rx^py, etc. 

Now, given that v does not depend on x and on y, the only non vanishing term in 
(1631) is the one containing dxdyW = B; so that we have 

k^yp = -dly wdlpV = -B dip V (64) 

o^y 3 

There is no contribution to the Ricci tensor coming from 9 Rp.xyp because the Minkowskian 

o^y 

metric tensor is diagonal and 9 =0; therefore, to make the Ricci tensor vanish up to 
the third order means that the following equation must be satisfied: 

9^'Kxyp = -B9^'dlpV = (65) 

This is the wave equation for v, imposing the dispersion relation a;^ = k"^. 

The particular form of the "modulating" function w{x,y) can be easily generalized. 
As a first step, let us consider an arbitrary quadratic form in variables x and y: let us 
put 

w = Ax'^ + B xy + C y'^ 
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3 

Besides the non vanishing component Rfj,xyp, two other components of Riemann tensor 
will appear now: 

3 3 

Rfixxp = 2^4, and R^yyp = 2C 
which have the same structure as the {x,y) component fl64|) : 

3 

R^xxp = -dl^ w dlpV = -2A dip v, 

kyyp = -dlyW dl^v = -2C dl^v, (66) 

The components (xx), {xy) and [yy) of the Ricci tensor vanish if the same condition 
fl65|) is satisfied; but now we shall also make sure that all other components of the Ricci 
tensor vanish, too, which will be true if the following trace is zero: 

9''''hxxp + hyR^.yyp = -dl^W dlpV-dl^W d^dpV 

= -{2A + 2C)dlpV (67) 

leading to the extra condition on the coefficients A and C, namely, A = —C, thus leaving 
only two degrees of freedom for the function w. This suggests the quadrupolar charac- 
ter of the gravitational wave, which deforms the space simultaneously in two directions 
perpendicular to the direction of propagation; notice that if w depended only on one 
transversal variable, say x, the vanishing of the Ricci tensor would impose w = (or a 
constant, which would not have any physical meaning at all). It is also worthwhile to 
note that the fact the planar wave solutions appear only at the third order of deforma- 
tion echoes the well known result obtained via linearization of the metric tensor, telling 
that gravitational waves are emitted when the third time derivative of the quadrupolar 
moment is different from zero. 

The same is true for any homogeneous polynomial of two variables x and y, provided 
it satisfies the two-dimensional Laplace equation d"^^ w + dyy w = 0. Finally, we can gen- 
eralize our result by stating that the deformation of Minkowskian space-time embedded 
as a hyperplane in an five- dimensional Euclidean ambient space leads to the vanishing of 
the Ricci tensor up to the third order in small parameter e if it has the form 

z^ = ee' + w{x, y) + 0(6^) (68) 

Provided that v satisfies uo"^ = c^k"^ and w satisfies the two-dimensional Laplace equation 
V^w = 0. 

3 

Taking into account that the corresponding Riemann tensor e R^^xp is linear both 
in V and w, we can compose by superposition a transversally polarized plane wave of 
arbitrary shape and spectrum, propagating with the phase velocity equal to the speed of 
light. 

The particular form of plane wave solution suggests also the form of a spherical 
wave. The first-order deformation far from the source should contain a factor propagating 
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in radial direction, while the second-order deformation should depend on the angular 
variables. We should not expect total vanishing of the second-order correction to the 
Riemann tensor like it happened in the case of plane waves. It is important that there 
will be no propagating terms at that order of approximation; static terms vanishing at 
spatial infinity like or can be neglected and in fact describe the approximation 
to the static part of the space-time deformation inevitably produced by the source of 
spherical gravitational waves. 

Let us start with the first-order deformation of Minkowskian space-time embedded 
as a hyperplane in some pseudo-Euclidean space; it has one component along one extra 
dimension perpendicular to the Minkowskian hyperplane M4. We suppose that is depends 
on the variables r and t only: 

y5 = (t, r) (69) 

Being perpendicular to the embedded manifold M4 as seen from the host space, this 

deformation does not contribute to the first-order correction to the Riemann tensor. In 

2 

order to evaluate the second-order correction to the Riemann tensor, i^^i^Ap, we need to 
insert the expressions for second covariant derivatives of f . In a fiat space parameterized 
by spherical coordinates the non-vanishing Christoffel symbols are: 

— r sin^ 9 
cote. 

and in the case when f is a function only of t and r the non-vanishing combinations are: 

VtVrV = dlv, VtVtv = dlv 

V<^V(p V = r sin^ OdrV. 

We consider the same deformation of the fifth coordinate as before: 

= evix^") + w{x'') + h{x^) + ... 

The first order deformation of the Ricci tensor is still zero. In order that the the 
second order be zero we have to satisfy the eq. fl571) . Because we are looking for radiative 
solutions, we shall neglect all terms which decay at spatial infinity more rapidly than 
1/r, keeping only the radiative part. If we set 

pi{ujt—kr) 

vix^') = 

r 

the only non-vanishing components of radiative character, i.e. behaving at infinity 
like 1/r are: 

2 2 ^—2i{Lut—kr) 

R^^ = sin^ eRee ~ ik sin^ 6 {k^c^ - u^) (70) 



-r 

,-1 



r9 



r: 



rtfi 



— sin 6 cos 6 
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and they do vanish provided that lJ^ = k'^c^. The third order correction to the Riemann 
tensor is given by the equation fl6T]) : iFiom this we can easily calculate the third-order 
correction to the Ricci tensor. Note that the third order deformation w has to depend on 
the angles 6 and 0, but must also have the dimension of a length. This is why we choose 

w{x^) =r Q{e,(f)) (71) 

the only non vanishing components (in radiative approximation) are: 

3 3 ~i(ujt~kr) 

R^^ = sin^ 6 Ree ik (sin^ 6 dlQ + 2 sin'^ Q + 9jg + cos Q sin Q deQ) . (72) 

Considering the Laplace operator acting on a function fl7T]) . we see that vanishing of the 
laplacian of w = r Q{6, ip) coincides with the condition of vanishing of the two non-trivial 
components of the Ricci tensor (1721) : with respect to the 9 and variables: 

= — 5— (sin^ 9 dgW + 2r sin^ 9 drW + d'^w + sin 9 cos 9 dgw) 

sin 9 ^ 

= ^ (sin^ 9d'^oQ + 2 sin^ 9 Q + dlQ + sin 9 cos 9 deQ) (73) 

r sin 9 ^ 

So the third order correction at the Ricci tensor vanish if w{x^) satisfy the Laplace 
equation, and the analogy with the plane wave solution we have previously found is 
complete: 



= e— - + e^w{r,9,(f)) + 0{e^) 



i(uit—kr) 



with 



V^w = 



The non-radiative terms behaving at infinity like and r""^ may give a hint as to the 
modifications of Schwarzschild metric that have to be made in order to compensate them 
thus solving the third-order Einstein equations exactly, or at least up to that order in the 
development in powers of r. This will probably suppose the existence of time-dependent 
variations of Schwarzschild background that would serve as the source of our spherical 
gravitational wave. 



4.2 Isometric embedding of Schwarzschild 's solution 

Isometric embeddings of Einstein spaces in pseudo- Euclidean flat spaces of various dimen- 
sions and signatures can be found in J. Rosen's paper in [12]. An embedding of the exte- 
rior Schwarzschild solution which is of particular interest to us , cited in Rosen's paper. 
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has been found by Kasner [I3|, who also proved that the embedding of Schwarzschild's 
solution in a five-dimensional pseudo-Euclidean space is impossible. Kasner's embedding 
uses a pseudo- Euclidean space with signature (+ H ) and is defined as follows: 



(- 


2MG\ 






cos 1 




r 








1 




2MG\ 




(- 




sin ^ 




r 





ct 



MG 
r 



2MG 



(74) 



z 



r sin 6 cos 
r sin 6 sin 
r cos 9 



Here M is the mass of the central gravitating body and G denotes Newton's gravita- 
tional constant. (Note the dimensional factor MG in front of the definitions of z^ and z"^ 
in order to give these coordinates the dimension of length) . 

The embedded four- dimensional manifold V4 is parameterized by the coordinates x^, 
with /i = 0, 1, 2, 3 so that 



X = ct, X 

Let us denote the flat metric by 
r]AB = diag (+ H 



r, X 



e, x^ 



(75) 



with 



A 5, 



1,2,. ..6. 



Then it is easy to check that the induced metric on the embedded manifold has indeed 
the usual Schwarzschild form: 



ds"^ = riAB d^z^ d^z^ dx'^dx" = g^^ (x^) dx^'dx" 

dr"^ 



(76) 



^ 2MG\ 2 , 2 

1 c^dt^ 

r J 



1 



2MG 



-r\de'-sm' ed^) 



However, this particular embedding is not unique. In 1959 C. Fronsdal [T^ proposed 
a similar embedding into pseudo-Euclidean space with the signature (l-|-,5— ), using 
hyperbolic functions instead the trigonometric ones. Fronsdal's embedding is defined as 
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follows: 



MG\ 1 



MG\ 1 



rl 



2MG\ 
r j 

2MG\ 
r j 

MG 



sinh 



cosh 



ct \ 

aJg) 

ct 
MG 



2MG 



dr, 



1 - 

rsin 6 cos 
rsin 6 sin 
r cos 6 



Again, it is easy to check that with the pseudo-Euclidean metric tjab 
the induced metric 



(77) 



diag{- 



g,j,u = Vab d^z^ djyZ 



B 



is the same as in (1761) . 

We have found a way to encode the two cases in a single formula. Introducing two 
constants a and x "we can write the pseudo-Euclidean six dimensional metric as 77^5 = 
diag{+a ), and the embedding functions z^{x^) as follows: 



2 exp 



MG\ 1 



MG\1- 



2MG 



2MG 



MG' 



2 exp 



X +X exp 



2X 



g-M4G4_2MGr3 
(2MG-r)r3 



(79) 



z'^ = r sin 9 cos 0, 
^5 _ j-gin ^ gin 0^ 

2-6 _ ^ cos ^ 

from which it follows that choosing the values (a = — l,x = 1) o^ig g^ts Fronsdal's 
embedding, while by choosing the values (a = 1, x = i) one gets Kasner's embedding. 
However in what follow we will use Fronsdal's functions. 

Let us also write down the non vanishing components of the Christoffel symbols and 
of the Riemann tensor, which do not depend on the embedding being inherent to the 
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internal geometry of Schwarzschild's solution. The Christoffel symbols are the following: 



_ _ MG ( 2MG\'^ 



h - -^(l-^)"' (80) 

MG 2MGy 



2MG\ 

-r 1 - 



r / 

= -r sm^^ 1 



On O 1 

ro ^ • 

er ~ J- T-6* ~ : 

r 



= -sin^cos^; 



) 



1 1 1 
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° ,n °w COS 9 , , 

= rL = — ^- 81 



~ sin^' 

while the non vanishing components of Riemann's tensor of Schwarzschild's metric 
are given by the following expressions: 



R 



o 



trtr 



2MG 



Rtete — 1 



MG f 2MG\ 



MGf^ 2MG\ . 2, 
Rt^tip = — ~ j^^"^ ^ 

MG r 2MGY^ , , 

Rr9r9 = U (82) 



MG f 2MG\~^ 



Rrprp = ~. j ^ 

Re^ep = -2MGrsm^9 

4.3 Deformations of the embedded Schwarzschild manifold 

As explained above, the transversality condition ( 12^ which represents four independent 
equations, shall leave only two arbitrary functions describing non-trivial deformations 
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of the exterior Schwarzschild solution in six-dimensional pseudo-Euclidean space. Let 
us find as simple choice as possible in order to make easier the subsequent calculus of 
Einstein equations. 

Let us examine the transversality conditions ( 12^ one by one in the concrete case of 
the embedding given by formulae ( 17^ . Only the first two embedding functions depend 
on time t; therefore the t-component of transversality condition becomes 



do + do 







1 



2MG\ 
r J 



sin 



ct 

Jig 



2 f 
V cos I - 



MG 



(83) 



It follows that the two functions and v"^ are proportional to a common function f (x^ 
and the ansatz that solves (l83l) is: 



= v{x^) z^ = MG vix^") sinh 
y2 = y(x^) = MG v(x^) cosh 



■'(—) 

\mg)' 

\MG J 



The common factor MG in front of the definitions above is put there to give the 
proper dimension (length) to and f ^. The unknown function u is then dimensionless. 
By analogy one can observe that the only embedding functions which depend on angular 
variables 6 and are the last three ones, z^, z^ and z^. Therefore the simplest way to 
satisfy simultaneously the last two equations corresponding to the components 6 and 
of (El) is to set 



^5) 



or more explicitly, 
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w{x'^) 

MG . ^ 
r 

MG , 
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= MG 
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v' = 




= MG 
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= MG 
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cos 


e. 



Again, we have conserved the functions z"^ in these formulae in order to maintain the 
spatial dimension in the definition of embedding functions and their deformations, while 
the unknown functions v and w will be kept dimensionless. 

Now we have two independent functions of four- dimensional coordinates, u{x^) and 
w{x^). The only remaining component of which we shall denote by = (MG) h{x^), 
is entirely determined by the radial component of the transversality condition fl24|) . 
Vab v^dr z^ = 0. It reads: 



M'^G'^vix") + hix^'y 



1- 



AfG\4 



2MG' 



+ w{x'') = 0. 
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After simplifying by the common factor {MG) we can express the function h[x'^) as the 
following linear combination of v and w: 



h{r) 



(M- 




So finally we can write 




(86) 



L r3 r-2MG )\ 



V 



,6 



MG sin 9 cos 0, 

MG w(x^) sin 9 sin 0, 
MG w{x^) cos ^. 



4.4 First order approximation of Einstein equations 

The system (15T!) represents an approximate version of the full system of Einstein equations 
with a source corresponding to the presence of extra matter in the vicinity of the central 
spherically symmetric mass. The extra matter gives rise to the energy-momentum tensor 
which can be that of of a test particle turning around the central body, or to an axially- 
symmetric distribution of matter, e.g. a homogeneous disc turning around the central 
body in the equatorial plane. It can be also a spherically symmetric electromagnetic 
field defined everywhere in the space surrounding the central body. The last solution to 
the problem is well know, it is the Reissner-Nordstr0m metric generated by a spherically 
symmetric mass endowed with an electric charge Q. 

In all these cases we should suppose that the absolute value of the energy-momentum 
tensor is very small when compared to the central mass M; in other words, the total mass 
m of the rotating disc or that of the point-particle representing a satellite should be small 
enough in order to treat the ration m/M as an infinitesimal parameter of the deformation 
of the background space-time provoked by its presence. All non-vanishing components of 
Riemann curvature tensor (and those of the Ricci tensor, too) in the case of Schwarzschild 
metric are proportional to M; in vacuo Schwarzshild metric is Ricci-flat and its Einstein 
tensor vanishes; but if the right-hand side is not zero anymore but proportional to m, 
Einstein's equations can be divided on both sides by M, and the right-hand side will 
become proportional to the dimensionless small parameter e = {m/M) 

If a solution of the linearized problem can be found, it will automatically generate 
second and higher order corrections on both sides. One should keep in mind the fact 



23 



that the energy-momentum tensor chosen as the right-hand side of Einstein's equations 
must be also a solution of the systam of equations ruling the motion of matter in a 
given background, e.g. the geodesic equation when one wants to describe the motion of 
a satellite with a very small mass in a given space-time geometry. But if this geometry 
changes (in our case via a perturbation proportional to the small parameter e), so should 
also the geodesic equations, because the connection coefficients have been modified by 
addition of small perturbation proportional to e. The geodesic equation thus deformed 
will lead to slightly modified solutions, the unperturbed one plus the perturbation linear 
in e. The source changed in this way, the right-hand side of Einstein's equations would 
contain now the terms proportional to m/M and the new ones proportional to e [m/M) = 
{m/My, and the second-order approximation of Einstein's equations should be solved, 
too, and so forth. 

Before trying to find a solution with a non vanishing right-hand side we should de- 
termine the solutions of the homogeneous system fB7|) . There is one solution to the 
problem which is quite obvious: the functions (x^) corresponding to an infinitesimal 
variation of the mass parameter, M — > M + SM. The new embedding defines just a new 
Schwarzschild solution with a slightly different mass. We have 

z^{x^,M + 6M)=z^{x,M) + ^-^^ + ...^ 



where we have defined 



6M 
M 



and 



(^') = ^dM- ^^^^ 
By differentiating the embedding functions with respect to the mass parameter M we 



24 



find (179 



1 _ 3MG . I 7777 . 

(x^) = MG ' = sinh ^ - ^/l - 2—^ ct cosh 

'l - 2 " ' 



MG MG S r MG 

r 



(a;^ = MG^^^L cosh ^ - . /l - 2^ ct sinh 



^_2MG MG V r MG 

r 

t;'^ (a;^) = / , = dr (88) 

V r3(-2MG + r)2y^gf 

= 
v^ixf") = 
= 

In order to find other solutions we have to consider the ten equations ( H71) . We search 
for f (x^) and w(a;^) on the form 

vit,r,e,(l)) = fM9v{r)K{4>)K{t) 

w(t,r,9, 0) = fw{0)gw{r)h^{(f))k^{t). 
1 

Consider first the 9(f) component of R 

jlg^ = (-M^GV -2r5)-ix 

X (^2rcot9{-r'^ + M^G^r~3M*G^)d^v 

+MG cot 6 [-r^ - M'^G^{3MG - 2r)] (9^w 
-2 r(-r4 + M^G^r - 3M^G'^)d^gV 

-MG [-r^ - M^G^{3MG - 2r)] OIqW 
this expression will vanish automatically if we set 

Me) = UO) = sme (89) 
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With this choice we get the following expressions for the t9 and 1 components 
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They are both equal to zero if we set 

kv{t) — ky](t) 

hv{<P) = K{(t)) 

y-vy ) ~ 3M^G^(2MG-r)ra ) 

The components r 9 and r (j) can be now written as 
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M*G^ -2MGr + r^ 



{3MG - r)g^{r) + {2MG - r) rg'^{r] 



which are both equal to zero if we set 
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l_2MG 



gUr) = (95) 



r 



Finally, the diagonal components are all zero if for the remaining functions we set: 

hw{<f>) = acos0 + 6sin0 

So finally we can write the deformation functions as: 



r^-SM^G^ + rM^^G^ L 2MG . ^ 

= 3r^M^GH2MG-r) V^-^""^" 

ct ct 
X (a sin (h + b cos 0) ( A sinh — + B cosh ■ 



MG MG 



^ _ 2MG 



w{x^) = MG- sin^ (asin0 + 6cos0) X 

r 

However, there is an essential difference between these two solutions. While the first 
one describes Riemann tensor whose all non vanishing components tend to zero at space 
infinity as r~^, the second one, although strictly Ricci-flat, has some components of 
Riemann tensor which do not vanish at space infinity; in the case when the hyperbolic 
functions are chosen we have: 
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r->oo MG 

Such a situation occurs in electrodynamics, too, where the constant electric or mag- 
netic field tensor obviously satisfies Maxwell's equations, but can not be considered as a 
regular solution at infinity. Nevertheless such field can be used as an independent solu- 
tion together with a regular solution, e.g. the Coulomb-like field, in a bounded portion 
of space. 
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In our case it is clear that the only solution acceptable at space infinity is the first 
one, generated by infinitesimal mass variation. This is similar to what occurs in classical 
electromagnetism: at very great distances all static fields generated by finite distribution 
of charges look at the first approximation as a Coulomb potential of the total charge; 
then comes the potential of a dipole momentum, etc. 

Having this second independent solution by no means could span the space of all 
solutions of our linear system. In the case of partial differential equations the space of 
solutions is infinite-dimensional, and in order to single out a solution one has to specify 
the initial or boundary conditions on an entire three-dimensional submanifold. 

Nevertheless the existence and the character of the two independent solutions of 
Einsteinian deformations of the Schwarzschild embedding, we conclude that any small 
perturbation of Schwarzschild solution, confined in a finite space surrounding the central 
body, is seen at infinity just as a Schwarzschild metric corresponding to the sum of the 
two masses, M + m. 

In a four dimensional space-time the above two solutions can coexist only separated by 
a three dimensional manifold dividing the four dimensional manifold in two disconnected 
parts. Such a discontinuity cannot correspond to a point-like object moving in the vicinity 



of the central mass along the worldline given by (s), because the Dirac delta function 
describing such an object with the following energy-momentum tensor: 

r^^*^ = m 5^{^''- {s)j M^(a;") u''{x'') (97) 
where the four-dimensional delta function means literally the tensor product: 



d^i^""- x'' (s) J = 6{t- t {s)) 6{r- r {s)) S{e- 9 (s)) 5(0- (s)), 

with u^{x") denoting the components of the four-velocity of the mass m. 

In this case the discontinuity is concentrated on a one-dimensional line which cannot 
divide the four dimensional space-time in two disconnected parts. 

However, the two solutions can be separated by e.g. a spherical surface which in 
the space-time is multiplied by the time axis, thus separating the space-time in two 
disconnected parts by means of Heaviside's function: 

till (^«) 

It is clear that far from the source only the well behaving solution survives, corre- 
sponding just to the variation of central mass without changing the symmetry. Other 
deformations can be found only if a definite symmetry breaking is chosen. 
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5 Second order deformations and axial symmetry 



5.1 General considerations 

The most important difference between the first and all the subsequent orders of defor- 
mation of embedded manifolds is that one can no more neglect the possibility of entering 
extra euclidean dimensions. The initial embedding in Ej^ can always be considered as a 

hyperplane in some higher-dimensional Euclidean space, with p' + q' — p — q = n. 

As it was shown in Section 3.1., the first order deformations going beyond this hyperplane 
do not contribute to the modifications of Einstein's equations. But starting from the sec- 
ond order in infinitesimal parameter e the deformations along the initial hyperplane and 
those from the extra dimensions will contribute simultaneously to the second-order mod- 
ifications to the Einstein equations, and the number of unknown functions to be found 
will become higher than before. 

The Kerr- Newman metric in Schwarzschild-like coordinates (see |21]) can be written 
as follows: 

2 , 2 AaMGr sin ^ , , S , 2 
: at H catdio dr 



^ 2 A 2MGr\ 



-S..^-(A + ^^^^^|!±^)smW (99) 

where 

E = + cos^^, A = - 2 MGr + a^. (100) 

It tends to the square of the Schwarzschild line element when a 0. 

It is natural to ask whether this metric can be also embedded E^^ , as was the case 
of the exterior Schwarzschild solution. One cannot exclude a priori such a possibility; 
the Reissner-Nordstr0m metric 



2^ ^-l 

0' 



2MG Q^G 1 2 2 7/i2 2 • 2 
h ^ I dr'' - r^d0^ - r^ sm"' 



can be embedded in Ef^ in the same manner as the Schwarzschild solution, but this 
is due to its spherical symmetry. 
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5.2 The Kerr metric as a deformation of Schwarzschild back- 
ground 

The axially symmetric Kerr-Newman solution can not be embedded in six (pseudo) Euclidean 
dimensions; it requires more dimensions than that, for the following reason. Whatever 
the dimension of ambient flat space, the embedding functions for this metric must de- 
pend on two extra parameters, ? and a. For each value of these two parameters we 
have a non-trivial four-dimensional Riemannian manifold embedded in a fiat space; but 
if we vary M and a independently, we get a two-parameter congruence of embedded 
four-dimensional manifolds, which can be looked upon as a six-dimensional manifold, 
which is by no means fiat. This can be checlced by calculus of the Riemann tensor of 
the corresponding six-dimensional metric. Therefore a higher than six-dimensional fiat 
space is required for the embedding of this non trivial six-dimensional manifold. 

In the early eighties R.R. Kuzeev has defined an isometric embedding of Kerr's metric 
( |99|) in a nine dimensional pseudo- Euclidean space ([22]). 

Before finding the embedding, it is much easier to define four Pfafiian one-forms whose 
squares summed up with appropriate signs will give the desired line element in diagonal 
form: 

ds^ = laJq — ujI — ujI — 0^1, (101) 

with 

1 

2MGrV , 2MGrasiY?e ( 2MGr , 
ujq = 1 — cat H — 1 — — dip, 



2 



^rfr, uj2 = V^de, ^3 = VAI 1- ) sine dr. (102) 



V 



with fllOOl) defined as above. 

It is easy to see that the squares of these differential expressions yield the Kerr- 
Newman metric when inserted in the formula fllOll) . and give the Schwarzschild line 
element in the limit a ^ 0. 

The embedding problem can be formulated now in terms of the integrability of these 
forms. One should find AN functions Q^'^{x^) satisfying the following identities: 

--dx^ = Q^^{x')uj, (103) 

The result is as follows: 
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2 I 2 
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(MG) I sin^sinf<^- 
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2ar 



smy cos I (fi 



MG 

ct 
MG 



(104) 



with functions $i (r, ^) and $2 (?", ^) defined by quite complicated integrals. 

When the angular momentum parameter a tends to zero, we recover the Schwarzschild 
metric embedded in a six-dimensional pseudo-Euclidean space, the last three embedding 
functions disappearing with a = 0. Yowever, the limit is not analytic because of the pres- 
ence of the square root of a in the formulas fllOSp . and radically different topology of the 
two solutions. Nevertheless we believe that this embedding may serve as a starting point 
for deformations describing non axially symmetric departures from the Kerr-Newman 
metric. This will be the subject of subsequent papers. 



6 Discussion and conclusions 

In this article we have set forth a new formalism that makes it easier to consider small 
perturbations of a given Einsteinian background without unphysical degrees of freedom 
that mar traditional computations based on the deformations of the metric. Here the 
embedding provides us with clear geometric criterion selecting physical degrees of freedom 
and eliminating the unphysical ones. 
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We have succeeded the construction of wave solutions in flat space, or in an asymp- 
totically flat Schwarzschild manifold at spatial infinity. These solutions can display any 
imposed form due to the possibility of linear superposition of Legendre polynomials. Once 
a radiative solution is chosen, we can extrapolate it towards the smaller values of r where 
the non-radiative terms prevail. These can be seen as the corrections to Schwarzschild 
metric close to the central body that are responsible for the emission of gravitational 
waves detected at spatial infinity. 

The treatment of this problem is the subject of the work in progress. 
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